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Abstract. 

Some alternative gravity theories allow the Universal matter distribution to 
single out the existence of a preferred frame, which breaks the symmetry of local 
O ■ Lorentz invariance (LLI) for the gravitational interaction. In the post-Newtonian 

paramctrization of semi-conservative gravity theories, LLI violation is characterized 
£h ■ by two parameters, a\ and «2- In binary pulsars the isotropic violation of Lorentz 

invariance in the gravitational sector should lead to characteristic preferred frame 
effects (PFEs) in the orbital dynamics, if the barycenter of the binary is moving relative 
to the preferred frame with a velocity w. 

For small-eccentricity binaries, the effects induced by 6l\ and a.2 (the hat indicates 
possible modifications by strong-field effects) decouple, and can therefore be tested 
independently We use recent timing results of two compact pulsar- white dwarf binaries 
with known 3-dimensional velocity, PSRs J1012+5307 and J1738+0333, to constrain 
PFEs for strongly self-gravitating bodies, by assuming the isotropic cosmic microwave 
background to single out a preferred frame. The time derivative of the projected semi- 
major axis is used to constrain a precession of the orbital plane around w due to PFEs. 
From this we derive a limit \6t2\ < 1.8 x 10~ 4 at 95% confidence level, which is the most 
constraining limit for strongly self-gravitating systems up to now, however still three 
orders of magnitude weaker than the best Solar system limit for the corresponding 
weak- field parameter a-i. 

Concerning &i, we propose a new, robust method to constrain this parameter, which 
avoids the probabilistic considerations inherent in previous methods. This method is 
based on the fact that a PFE-induced intrinsic eccentricity cannot stay unobserved 
during a long-term observation due to the significant precession of periastron in binary 
pulsar with short orbital periods. Our most conservative result, a.\ = — 0.4jl|' ^ x 10~ 5 
at 95% confidence level from PSR J1738+0333, constitutes a significant improvement 
compared to current most stringent limits obtained both in Solar system and binary 
pulsar tests. 

We also derive corresponding limits for 6l\ and ui for a preferred frame that is at 
rest with respect to our Galaxy, and preferred frames that locally co-move with the 
rotation of our Galaxy. 

These limits will continue to improve significantly with future pulsar timing 
observations conducted at large radio telescopes. 



PACS numbers: 04.80. Cc, 11.30. Cp, 97.60.Gb 
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1. Introduction 

Lorentz invariance is one of the most important ingredients inherent in modern 
theoretical physics, including the standard model of particle physics and general 
relativity (GR). From a group theoretical viewpoint, it composes of two parts, rotational 
invariance and boost invariance. Rotational invariance forms a compact group, i.e., the 
SO(3) group, which can be probed throughout, while boost invariance forms a non- 
compact group, hence it, in principle, cannot be tested thoroughly, and deserves more 
scrutinies. Lorentz invariance is examined in particle physics to high precision [30], while 
not so well tested in gravitational physics [53], due to the challanges in gravitational 
precision experiments. 

On the other hand, some alternative gravitational theories predict the existence 
of a preferred frame, which might be singled out by the matter distribution in our 
Universe or through historical relics of vectorial or tensorial vacuum expectation values, 
if gravitational interaction is mediated by a vector field or a second tensor field, in 
addition to the canonical second-rank symmetric tensor field [51] . These theories include 
vector- metric theories |54[ [52] . Einstein- iEther theories [22], TeVeS theories [5], and 
standard model extensions of gravity [3]. 

In the parametrized post-Newtonian (PPN) formalism, preferred frame effects 
(PFEs) are characterized by three parameters, ax, «2 and «3 [HH [52]. Since «3 also 
causes energy-momentum conservation violation, and is well constrained to very high 
precision (see e.g. [48J, where the strong-field version of \as\ is constrained to be less 
than 4 x 10~ 20 at 95% confidence level, by using wide-orbit binary millisecond pulsars), 
we will not consider it further in this work. Experimental tests on PFEs induced by ax 
and «2 are roughly divided into three catalogues, i.e., geophysical tests, Solar system 
tests, and pulsar timing tests. 

Nordtvedt and Will [13] derived possible experimental indications of a preferred 
frame for the gravitational interaction in geophysics and orbital motions, e.g., an 
anomalous 12-hour sidereal tide of the solid Earth, an anomalous yearly variation in 
the rotational frequency of the Earth, and an anomalous perihelion shift of the planets. 
By now, the best limit for a\ in the Solar system comes from Lunar Laser Ranging 
(LLR), that places a 95% confidence limit of a\ = (— 0.7 ± 1.8) x 10~ 4 [39]. For the a2 
parameter, Nordtvedt [12] used the close alignment of the spin axis of the Sun and the 
total angular momentum vector of the Solar system to limit jo^l < 2.4 x 10~ 7 , under 
the assumption that the above two vectors were aligned when the Solar system formed 
five billion years ago (note, a Nordtvedt _ l a<2 y Damour and Esposito-Farese developed 
a method to put tight constraints onto the strong- field counterpart of ax, namely &x, 
from timing experiments of small-eccentricity binary pulsars [12]. Their calculation 
shows that the observational eccentricity vector, e(t), is a vectorial superposition of a 
"rotating eccentricity" e#(t) with constant length e#, and a fixed "forced eccentricity" 
ep. From probabilistic consideration, they were able to constrain |di| to be less than 
5.0 x 10 -4 (90% C.L.). This limit has been improved by a factor of three in [6], based on 
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the small-eccentricity binary pulsar PSR J2317+1439. Wex [5U] extended this method 
in statistically combining multiple systems, by taking care of a potential selection effect 
when simply picking the system with the most favorable parameter combination. He 
got a slightly improved result of \ai\ < 1.2 x 10~ 4 (95% C.L.). In section HJ we extend 
the statistically dependent method into a robust one which not only avoids involving 
probabilistic considerations concerning certain unobservable angles, but also gets a new 
constraint, a.\ = — 0.4jlQ x 10~ 5 (95% C.L.), that surpasses the current best constraints 
of both weak and strong fields. 

For the a<i parameter, because of the tight limit of [12], Damour and Esposito- 
Farese p2j dropped the a 2 term when calculating binary orbital dynamics. We 
stress that the limit of [12] is obtained in a weak-field gravitational environment, 
while in the strong- field regime, like inside a neutron star (NS), a 2 might take 
an independent value, significantly different from the Solar system value. In fact, 
it was discovered that in certain classes of tensor-scalar theories of gravity, large 
non-perturbative strong-field deviations from GR can occur, through a phenomenon 
called "spontaneous scalarization" [U]. Although tensor-scalar theories of gravity 
are conservative gravity theories and do not show PFEs, it is natural to assume the 
possibility of similar non-perturbative effects in gravitational theories with local Lorentz 
invariance (LLI) violation. Hence, we feel that it is still worth to independently test, in 
a phenomeno logical approach, the strong-field counterpart of a 2 , namely a 2 , in pulsar 
binary timing experiments. 

Following [12], we calculate the a 2 effect for pulsar binaries, and find that it 
practically decouples from the ati effect for small orbital eccentricities (e <C 1). The 
otx term tends to polarize the eccentricity vector towards a direction perpendicular to 
the orbital angular momentum and the binary barycentric velocity with respect to the 
preferred frame, w [12]. It causes dynamical effects inside the orbital plane. In contrast, 
the a 2 term imposes a precession of the orbital angular momentum around the direction 
of w. It causes a change in the orbital inclination angle with respect to the line of sight, 
i. Consequently, a 2 induces a non- vanishing time derivative of the projected semi- major 
axis. 

Wex and Kramer [51] developed a pulsar timing model that includes PFEs, by 
extending the widely used D amour- Deruelle timing model [TP] . Based on this model, 
they analyzed the time of arrivals (TOAs) of the double pulsar, PSR J0737— 3039A/B 
[H E3 EU [321. and jointly limited &i and a 2 to be, —0.5 < a.\ < 0.3 and —0.3 < a 2 < 0.2, 
respectively^ Their analysis utilized two specific aspects of the double pulsar: 1) The 
measurement of the mass-ratio via the "double-line" nature of the system, and the 
measurement of the Shapiro delay allowed for a theory-independent determination of the 
effective gravitating masses of the two pulsars; 2) The large rate of periastron advance, 
Cj = 16.9degyr _1 [SI], would significantly change the binary orientation with respect to 

I They use a\ and 0% instead of &i and &2 in their notation. The limit quoted above assumes the 
preferred frame to be at rest with respect to the isotropic cosmic microwave background (see the original 
paper for constraints on other directions). 
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a preferred frame within just a few years, leading to distinct long-term periodic effects 
in the orbital parameters. Consequently, the double pulsar has even the potential to 
measure di and a 2 , if they are non-zero. As emphasized in [51], double NS systems probe 
different aspects of a violation of LLI in the gravitational sector (interaction between 
two strongly self-gravitating bodies) from other kinds of binaries, e.g., NS- white dwarf 
(WD) binaries. Furthermore, from the simulations in [51] one expects that by now the 
precision of the PFE test with the double pulsar has greatly improved compared to the 
numbers in [51] . 

To this point, let us briefly summarize current best limits on the LLI violation in 
the weak field and strong field. 

(i) Weak field 

• From LLR [39], 

ai = (-0.7 ± 1.8) x 10" 4 (95%C.L.). (1) 

• From the close alignment of the spin of the Sun with the total angular 
momentum of the Solar system [32] , 

\a 2 \ < 2.4 x 10~ 7 . (2) 

One should be aware that this result depends on an assumption about the 
alignment of the spin of the Sun and the angular momentum of the Solar 
system right after their formation, five billion years ago. LLR experiments [39] 
get a weaker limit, a 2 = (1.8 ± 5.0) x 10~ 5 (95% C.L.). 

(ii) Strong field 

• From the population of small-eccentricity NS- WD binaries [50] , 

< 1.2 x 10~ 4 (95%C.L.). (3) 

• From a NS-NS system, namely the double pulsar [51] . 

- 0.3 <a 2 < 0.2 (95% C.L.) . (4) 

In this paper, we derive the full secular dynamical evolution of a pulsar binary 
system of arbitrary eccentricity, under the influence of both &i and a 2 . Afterwards, 
we utilize our analytical results to propose two new methods to constrain &% and a 2 , 
respectively. By using Monte Carlo simulations, we are able to get stringent limits from 
small-eccentricity NS-WD binaries, PSRs J1012+5307 [331 [M] and J1738+0333 [U|20], 
with measurement errors properly accounted for. 

The paper is organized as follows. In section [2] we derive the orbital dynamics 
from a generic semi-conservative Lagrangian keeping both the &i and the a 2 terms, 
for arbitrarily eccentric orbits. We find that in the limit of a small eccentricity, 
these two parameters decouple — 6t\ affects the evolution of the eccentricity vector 
in the orbital plane, while a 2 controls the precession of the orbital angular momentum. 
We introduce the isotropic cosmic microwave background (CMB) frame as the most 
important preferred frame, for our subsequent numerical calculations in sections [3] and HI 
In section|3l \a 2 \ is derived to be less than 1.8 x 10~ 4 (95% C.L.) from timing experiments 
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of pulsar binaries PSRs J1012+5307 and J1738+0333. In section 0], we develop a new, 
robust method to constrain d 1; which overcomes the need of probabilistic considerations 
inherent in the former methods. The most conservative limit, &\ = — 0.4l 3 \ x 10 -5 
(95% C.L.), is derived from the PSR J1738+0333 binary system. Section [5] gives the 
corresponding results on a\ and &2 when the Galaxy or the local Galactic rotation are 
assumed to single out a preferred frame. The limits for these frames of reference are 
found to be slightly weaker than the ones for the CMB frame. Furthermore, we give a 
discussion on the strong-field aspects of our tests, discuss future improvements of these 
two tests, and briefly summarize the results of the paper. 



2. Binary dynamics of the semi-conservative Lagrangian 



We consider the two-body dynamics of a binary system consisting of a pulsar with mass 
m p and its companion with mass m c . In the presence of a preferred reference frame, 
the orbital motion of such a system is described by a two-body non-boost-invariant 
Lagrangian [52} [12] 

L = Ls - + + L& 2 . (5) 

The Lagrangian fl5]) consists of Lg&, the post-Newtonian (PN) terms from GR and 
its minimal extensions characterized by the (strong-field) Eddington-Robertson-Schiff 
parameters, (3 and 7, 
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and the velocity-dependent, non-boost-invariant terms, related to non- vanishing &i and 

&2, 



Gm p m c 



(vj • v c °) 



2c 2 



«2 



Gm p m c 



:v°-v c )-(n. V 0)( n .v 6 ) 



(7) 

2c 2 ' (8) 

where M = m p + m c , r = |r| is the coordinate separation of two components, n = r/r, 
v° denotes the "absolute" velocity with respect to the preferred frame, and c is the 
speed of light. In the above Lagrangian, we add a "hat" onto the notations of 7, /3, 
ati, °-2i and also the gravitational constant G, to underline the fact that we are dealing 
with the PN parameters associated with compact objects, where strong-field effects 
might contribute to these values, making them different from their counterparts in the 
weak field. The specific dependence on the strong-field contributions depends on the 
gravitational theories under consideration. In GR, one finds G 

Oil = &2 = 0. 



G, (3 — 7 = 1, and 
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The Lagrangian (J5J) can be obtained from a more generalized, semi-conservative 
Einstein-Infeld-Hoffmann Lagrangiarjjj in the calculations of [52] and [51], by setting 



A p = A c = 1, g = G/G, B/G = |(27 + 1), C/g = |(4 7 + &! - a 2 + 3), £/£ = d 2 + 1, 
X c V p /g 2 + X p V c /g 2 = 2/3-1, where X p = m p /M and X c = m c /M. 



The assumption A p = A c = 1, made in this paper, requires additional justification. 
These parameters are equal to one at first post- Newtonian order [52], but for strongly 
self-gravitating bodies, e.g. NSs, they could significantly deviate from one (e.g. in 
Einstein-iEther theory [19]). The subject of this paper are pulsars with low- mass WD 
companions, and therefore A c — 1. Concerning A p , it enters the secular changes of the 
relevant orbital parameters only as addition to the parameters &i and a 2 (see [51] for 
details). In the first case the A p contribution is multiplied by a factor of 2X 2 and in 
the second case by a factor X c . Since X c ~ 0.1, for our binary systems, in both cases 
the A p contributions are expected to be small with respect to the c\\ and a 2 terms. 

Besides the PPN parameters (3, 7, a% and a 2 , semi- conservative theories of gravity 
could have a non-zero Whitehead term, characterized by £ [52]. It reflects preferred- 
location effects, such as an anisotropy in the local gravitational constant caused by an 
external gravitational potential. Various well-motivated (fully conservative and semi- 
conservative) gravity theories have £ = (for instance, see Table 3 of [53]). Therefore, 
we will only include (3, 7, di, and a 2 in our following discussion, and ignore a potential 
strong-field counterpart of the Whitehead term. We note in passing, that for small- 
eccentricity binaries, the presence of a Whitehead term only changes the a 2 test, which 
then constrains a combination of a 2 and £. The a\ test, on the other hand, remains 
unchanged. This can be readily seen from (8.73) in [52] . 

2.1. Orbital dynamics in the presence of PFEs 

Using the Euler-Lagrange equations and the post-Galilean transformations [9], we can 
get the relative acceleration for a pulsar binary system in the comoving frame, whose 
center of mass moves relative to the preferred frame with a velocity w, 



where A^ N ^ is the "Newtonian" acceleration, A^ N ^ = — GMn/r 2 , A( PN )/c 2 is the first 
PN acceleration without w-dependent contributions, and A^/c 2 is the additional 
acceleration from the motion of the binary system with respect to the preferred frame. 
For expressions of these accelerations, see [52l [12], [5T] • 

As is well known, the acceleration A^ PN ^/c 2 produces a secular advance of the 
longitude of periastron, 



§ A "semi-conservative" Lagrangian corresponds to a gravity theory that possesses conservation laws 
for the total energy and momentum. Any theory that is based on an invariant action principle is 
"semi-conservative" [S3] . 



r = r p _ r c = A W + A (pn) f(? + A (w) j(? ? 



(9) 




(10) 
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Figure 1. Illustration of the geometry of the binary system, and the notation used in 
the paper. The coordinate system (Iq,Jo,Ko) corresponds to (Iq, Jo, Kq) in figure 1 
or (Io,Jo,Ko) in the text of [18], and the coordinate system (a, b, k) corresponds to 
(a, b, c) in [T2]. w is the velocity of the binary system with respect to the preferred 
frame, while wi is its projection into the orbital plane. 



where 

Vo = {GMn b ) l '\ (11) 

T = ^(2 + 27-/3) + ^(2d 1 -a 2 )X p X c , (12) 
3 o 

and rib = 27r / 'P& is the orbital frequency of the binary systemijj The "characteristic" 
velocities in GR, V^ GR) , for pulsar binaries PSRs J1012+5307 and J1738+0333, are given 
in table[TJ In GR, J 7 = 1. Because A ( - PN ^/c 2 lies in the orbital plane, it has no effect on 
the longitude of the ascending node Q, and the orbital inclination angle % (see figure [T] 
for illustration of these angles). In addition, A^ PN ^/c 2 is verified to have no effect on e, 
the length of the eccentricity vector, and a, the semi-major axis of the relative orbit (cf. 
(115]) below, with w = 0) 

As for the acceleration from PFEs, A^/c 2 , Damour and Esposito-Farese [12] 
worked out the influence of cii-related terms on the orbital evolution. After averaging 
over one orbital period, they found for changes in the semi- major axis a, 1 = y/1 — e 2 k, 

The velocity Vo, which corresponds to /?oc in [18], should not be confused with vq = Vo/Vl — e 2 
in (17-21) of [ST]. 
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and the eccentricity vector e = e a 
da \ 



dt I 



= 0, (13) 

PN+ai 



§) 4ki n6V ° efe(bxw) ' (i4) 

• U I -\ ; TP f ft 2 , K e2 ^ fc 



dt/ 



PN+Q! 



ew PN b + ^ — \n b V F e ( \/l - e 2 u> a a + u> b b 6 k ) , (15) 

where g = m p /m c is an observable quantity, due to the additional optical information 
for the two pulsar binaries used in this paper. The three unit vectors (a, b, k) form a 
right-handed triad of a coordinate system with its origin at the center of mass of the 
binary system, where a points to the position of periastron, k points along the orbital 
angular momentum, and b illustrated in figure [TJ Furthermore, 

F - S urjr=?- (16) 

is a function that, for bound orbits (0 < e < 1), takes a value in the interval [|, l). 

In addition to the contribution from the d\ and PN terms, we find that for the 
change caused by the a.2 terms, after averaging over one orbital period, 
da ' 
dt 



0, (17) 



dl \ 6>2 



■ 



dt/, c 



—n b F e (w k k + e 2 F e w b h) x w 



(^) = ^ n b p e (FeVl - e 2 w a w b a - F e Wa 2 Wb b + w b w k kj e , (19) 

where (w a ,Wb,Wk) are the coordinate components of w in the (a, b, k) system. 

From (|T5|) and (IT7|) we can see that, as expected, to first order there is no change 
in the semi-major axis of the orbit from PFEs. 

2.2. Small- eccentricity orbits and PFEs 

The coupled differential equations above simplify considerably for small eccentricities. 
When e C 1, one finds F e ~ 1/2 and 1 ~ k. To leading order in the (numerically) 
relevant contributions, (I14p and (ITSj) become 

f) =">• < 20 > 

dk\ d.2 



dt/, - 2c2 n b w k kxw, (21) 
and f|T5|) and ffT9]) simplify to 

— ) ~e WpN b + — — — n fe V w ± , (22) 

ut / PN+&1 ^ L y ~r 1 

'de\ . . 
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where w = |w| and wj_ = w a dL + wt,h is the projection of w into the orbital plane. 
Above four equations have been derived, under the consideration that e < 10 -6 , 
w 2 /c 2 ~ Vovj/c 2 ~ Vq/c 2 ~ 10~ 6 , for the NS-WD systems which are to be used in 
our calculations, i.e., PSRs J1012+5307 and J1738+0333 (see tabled]). 

From (1201 and (1221) . Damour and Esposito-Farese [12J worked out the eccentricity 
vector evolution under the influence of the PN and 6t\ terms, which turns out to be 
a superposition of a PN- induced precessing eccentricity e R (t), and a constant "forced 
eccentricity" e^, introduced by a±. 

In terms of a geometrical interpretation of the time evolution of the orbital 
eccentricity, the physical consequence of 6>\ was extensively studied p21 El ED]- The 
physical consequence of a 2 for small-eccentricity binary systems is readily derived 
from ( 12TT) . which shows that a non-zero a.2 causes a precession of the orbital angular 
momentum around the fixed direction w with an angular frequency, 



where tp is the angle between k and w (see figure [T] for an illustration of the orbital 
geometry and the orbital angular momentum precession). To leading order, this 
precession is purely determined by a 2 (see (120]) ). 

The precession (I24p induces a secular change of the projected semi- major axis of 
the pulsar orbit. The rate of change is given by 



where $ is the angle between wi and the direction of ascending node. In section El we 
will apply (|25|) to constrain d 2 from two relativistic small-eccentricity NS-WD binaries, 
namely, PSRs J1012+5307 and J1738+0333. 

2. 3. The preferred frame 

As the most natural preferred frame for our following calculations, we choose the 
frame determined by the isotropic CMB, like this is generally done in the literature 
on preferred-frame tests. To use other frames, the generalization is straightforward. As 
an example, in section [5] we also present limits on ct\ and &2 for which the Galaxy or 
the local Galactic rotation is assumed to determine the preferred frame. 

From the five-year Wilkinson Microwave Anisotropy Probe (WMAP) operations, a 
CMB dipole measurement of 3.355 ± 0.008 mK was obtained, which implies a peculiar 
velocity of the Solar system barycenter (SSB) with respect to the CMB frame of 
I v ssb-cmb| = 369.0 ± 0.9 kms -1 , in the direction of Galactic longitude and latitude 
(l,b) = (263.99° ± 0.14°, 48.26° ± 0.03°) [21]. Results from the seven-year WMAP 
observations are unchanged [21]. The binary velocity with respect to the preferred 
frame is w = v P sr_ssb + v ssb-cmb, where v P sr_ssb is the 3-dimensional (3D) motion 
of the pulsar binary system, with respect to the SSB. For PSRs J1012+5307 and 
J1738+0333, vssb-cmb can be derived from a combination of the distance and proper 




(24) 




(25) 
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motion measurements from radio timing, and the radial velocity obtained from spectral 
observations of the WD. 

3. New constraints on a 2 

To constrain &2 from binary pulsar observations, we start from (1251) . where the orbital 
frequency, n& = 27r/P&, and the projected semi-major axis, x, are observable Keplerian 
parameters, while the time derivative of x, x, belongs to the set of phenomenological 
post-Keplerian parameters [10J [18]. They are obtained with high precision from radio 
timing observations. In (1251) . we also need the inclination of the binary orbit with respect 
to the line of sight i. For the binaries of this paper, i can be determined (modulo the 
ambiguity of i — > 180° — i) from the mass function, leading to 

sinz = ^(g + l). (26) 
vo 

The companion mass m c is inferred from spectroscopic and photometric studies of the 
WD companion using well tested atmospheric model for such WDs [U [1] . The pulsar 
mass m p is determined from the mass ratio of the pulsar and its companion, q = m p /m c , 
which is inferred from the radial velocity and the orbital parameters of the binary 
system [H [T]. Unfortunately, the information is not sufficient to calculate smi, since 
Vo contains, besides the known total mass M, the effective gravitational constant G 
which is a priori unknown if one does not specify a given gravity theory. In principle, 
strong-field modification could lead to a significant deviation of G from the GR value, G. 
Such modifications, on the other hand, are expected to be accompanied by a significant 
amount of dipolar gravitational radiation (as an example, see [T2J for the case of tensor- 
scalar theories of gravity), which is neither the case in PSR J1012+5307 [Hj nor in 
PSR J1738+0333 [20] . Consequently, for the required precision in i, we can safely 
assume G ~ G in ( 1261) . 

In order to fully determine the orientation of the binary with respect to w (if) and $ 
in ( 125]) ). one also needs the longitude of the ascending node, Q, an angle which (in most 
cases) is not measurable from pulsar timing experiments. Consequently, in our a<i tests 
we will treat Q as a random variable uniformly distributed between 0° and 360°. This, 
however, will require probabilistic arguments in order to exclude those (small) ranges 
of Q where d 2 would practically be unconstrained. 

3.1. PSR J 10 12+5307 

PSR J1012+5307 is a small-eccentricity NS-WD binary system, with an orbital period 
of ~ 14.5 h. The pulsar was discovered in 1993 with the 76-m Lovell radio telescope at 
Jodrell Bank [41], and optical observations revealed its companion being a helium WD 
|36| . Callanan et al [8] measured a systemic radial velocity of 44 ± Skms^ 1 relative to 
the SSB, the mass ratio q = 10.5 ± 0.5, and the companion mass m c = 0.16 ± 0.02 M Q . 

Lange et al [33] used 4-year Effelsberg 100-m radio telescope timing data and 
7-year 76-m Lovell telescope timing data to derive a set of pulsar timing parameters, 
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Figure 2. Upper: Illustration of two important angles of the PSR J1012+5307 binary 
in our calculation, as a function of the unobservable longitude of the ascending node fl: 
tp, the angle between w and the orbital angular momentum, and the angle between 
wx and the direction of ascending node. Middle: Proper motion contribution to i, 
namely x FM (gray), and the residual i, namely i obs — i PM (blue). Lower: Derived ct2 
from the residual x. Different contours correspond to 68%, 95%, and 99.7% confidence 
levels. The left column is for the i < 90° branch, and the right column is for the 
i > 90° branch. 
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by utilizing the low eccentricity binary timing model ELL1. These timing results were 
used to put stringent limits on the emission of dipolar gravitational radiation by this 
NS-WD system. Most recently, Lazaridis et al [M] updated the timing parameters 
by using 15 years of observations from the European Pulsar Timing Array (EPTA) 
network, consisting of the four radio-telescopes Effelsberg (Germany), Jodrell Bank 
(UK), Westerbork (the Netherlands) and Nangay (France). 

The short orbital period and the measurement of its spatial systemic velocity make 
the PSR J1012+5307 system particularly interesting for studies of PFEs and tests of the 
corresponding parameters. To utilize (125]) . the full information of w and the orientation 
of the orbital plane is needed. To get w, we calculate the binary velocity with respect 
to SSB from the transverse velocity obtained from radio timing and the radial velocity 
measurement via spectroscopy of the WD companion. With these measurements at 
hand one can compute w = vpsr-ssb + vssb-cmb- The inclination of the orbital plane, 
i, is calculated from (|26l) . with a sign ambiguity between i < 90° and i > 90°. Hence 
we have two branches of solution. Moreover, pulsar timing experiments generally give 
no information on the longitude of the ascending node Q. We sample it in the range 
[0,360°). Figure [2] illustrates ip (angle between w and the orbital angular momentum) 
and $ (angle between the ascending node and w_|_), as a function of Q. In the figures of 
ijj and i?, measurement uncertainties of m c , q, and the proper motion are not included. 
However, our simulations to constrain 6:2, which are to be discussed below, take full 
account of all measurement uncertainties. 

To look into the change of x induced by the 612 term, we should separate other 
potential effects from the measured x. A change of x can come from various astrophysical 
and gravitational effects [18] [35]. All effects that cause a change to the semi-major 
axis of the system, like gravitational wave damping and component mass loss, can be 
constrained observationally via the observed P b obs (see tabled]). In fact, we can re- write 
(8.76) in [35] to 

The remaining terms are due to a change in the Doppler factor D, the proper motion 
of the binary system, a change in the aberration due to a change in the pulsar-spin 
orientation, a change in the orbital inclination due to spin-orbit coupling effects, and 
finally a mass distribution in the vicinity of the system. Before calculating 612 from ( 125]) . 
these influences on x should be subtracted. We will discuss them term by term in the 
following. 

Using the measured quantities of PSR J1012+5307 in table [1] the first term on the 
right-hand side of (12?"]) can be estimated to be x Pb ~4x 10 -19 s s _1 , which is four orders 
of magnitude smaller than the relevant scale. 

The second term, —D/3D, includes contributions from the Galactic acceleration of 
the binary system, and the Shklovskii effect that is induced by the transverse proper 
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motion 146 1. One finds 



D 1 v 2 
- — = -Ko • (gPSR - gss B ) + > (28) 
D c ca 

where K is the unit vector pointing from the SSB to the binary pulsar, as defined before, 

al is the pulsar distance from the SSB, = d^J ju| + ju| is the transverse velocity of the 

system with respect to the SSB, gpsR and gssB are the Galactic accelerations of the 

binary and the SSB, respectively. The contributions from the Galactic acceleration and 

the Shklovskii effect to x are of order 2 x lCT^ss -1 and 3 x 10~ 19 ss -1 , respectively. 

Hence, both contributions are negligible |34j . 

The third term of (|27|) is a variation of x caused by a change of the orbital inclination 

i, due to the proper motion of the binary system j2j [29] . 

' . \ PM 
X \ 

= (— j2 a sin f2 + fig cosf2) cot i . (29) 

x J 

This contribution is not negligible. The contribution from the proper motion effect 
is depicted in the middle panels of figure [2] for i < 90° (left) and i > 90° (right), 
respectively, as a function of the unknown longitude of the ascending node Q. In [31] 
the x measurement was used to constrain Q, by assuming that the measured x is solely 
caused by the proper motion effect ( J29l) . They got constraints on Q, by requiring 
x obs = x FM (see (11-14) in [31]). In our test we have to keep full ignorance of Q, and 
cannot assume x obs = x PM . The residual x, after subtracting of the contribution from 
binary proper motion, is also plotted in the same figure. The intersections of x obs — x PM 
and the horizontal null lines in the figure correspond to the limit on Q obtained in [34J. 

The fourth term in (1271) is due to the varying aberration caused by geodetic 
precession of the pulsar spin axis [18]. For a nearly circular orbit one finds 

&Sa P cot A sin 2r/ + cot i cos n . . 

~TT — FT ■ 7 "good' 30 

dt Pb sin A 

where P is the pulsar spin period, A and rj are positional angles of the pulsar spin vector 
(see figure 1 in [18] for details). In GR, the geodetic precession rate for a nearly circular 
orbit is given by [1] 

n 3 + 4g_ M GR) V 

n ^^2tTTW [~^J nb - (31) 

Consequently, from the timing parameters in table [T] one finds that the de^/dt term 
produces a change of x of order 10 _18 ss _1 for typical spin orientations. Hence, it is 
negligible in our case, unless there is a deviation from GR by at least a factor of 100, 
which we consider as highly unlikely, as such a large deviation of gravity in this system 
is clearly not seen in the gravitational wave emission [31]. Moreover, PSR J1012+5307 
is a highly recycled pulsar, and therefore its spin axis is expected to be nearly aligned 
with the orbital angular momentum (i.e. rj ~ —90° and A ~ i), which greatly suppresses 
this effect anyway. 
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The classical spin-orbital coupling, due to the quadrupole moment of the companion 
star, can change the inclination of the orbital plane, which induces the fifth term of ( 127|) . 
But this effect is only important for main-sequence star companions [261 H9] or rapidly 
rotating WD companions [27], and can be neglected here. 

The last term of (1271) is only valid if there is a nearby third companion, that 
perturbs the orbit significantly, like in the PSR B1620— 26 system [2j. This is not the 
case for PSR J1012+5307 [M|, as this would be well seen in the presence of higher oder 
derivatives of the rotational frequency of the pulsar [25] . 

In summary, only the proper motion term (129]) is important in our studies here. 
After accounting for this effect on x ohs for every given Q, we calculate the contour plots 
of a 2 , and present them in the lower panels of figure |2] for i < 90° (left) and i > 90° 
(right), respectively. In the calculation, 10 5 Monte Carlo simulations are implemented 
to account for the measurement uncertainties of [i a , [is, d, v r , q, m c , and x. As we 
can see, a 2 can be constrained to the order of 10~ 4 for most Q realization. The Q 
values in the figure where a 2 is virtually unconstrained correspond to the configurations 
when the angle between w and the orbital angular momentum if) ~ 90°, or the angle 
between the projected w onto the orbital plane wi and the direction of the ascending 
node d ~ 90° (see horizontal lines in the upper panels of figure [2J for corresponding 
angles and compare them with the divergencies in the lower panels). The reason for the 
divergencies is easy to see from ( 1251) . where the right hand side vanishes when ip = 90°, 
or $ = 90°, independent of a 2 . In this situation a 2 cannot be constrained. 

3.2. PSR J 1738+0333 

PSR J1738+0333 is a small-eccentricity NS-WD binary system, with an orbital period of 
~ 8.5 h, which, together with other well measured physical quantities, makes it a superb 
astrophysical laboratory to test gravitational theories [20]. The pulsar was discovered 
in 2001 in the Parkes high Galactic latitude survey [23], and later regularly timed with 
the 305-m Arecibo telescope [20]. It is one of the four millisecond pulsars known to be 
orbited by a WD companion bright enough for detailed spectroscopy [1] , among which, 
PSR J 1738+0333 is the most relativistic. Detailed optical studies of the WD companion 
and radio timing studies of the pulsar are presented in pQ and [20] , respectively. Thanks 
to their studies, accurate binary parameters and spatial motion (transverse and radial) 
are available for the PSR J1738+0333 binary system. Therefore, it also presents a good 
laboratory to study PFEs. 

The strategy to constrain a 2 is the same as in the case of PSR J1012+5307. First, 
we get w from vp SR _ SSB and vssb-cmb, and then the configuration of the system with 
respect to the CMB frame can be obtained, as a function of Q, with a sign ambiguity 
of i. Two important angles, ip and are depicted in the upper panels of figure El for 
% < 90° (left) and % > 90° (right), respectively. 

Along the line of arguments for PSR 1012+5307 one finds, using the results of [20] 
(see table [1] for binary parameters), that also for PSR J1738+0333 the only relevant 
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Figure 3. Same as figure El for PSR J1738+0333. 

term in ( )27|) is the contribution by the proper motion of the system, x PM . Because of 
the more precise measurements of the PSR J1738+0333 parameters, the uncertainty of 
x PM is smaller correspondingly, as illustrated in the middle panels of figure [3] for i < 90° 
(left) and i > 90° (right), respectively. The corresponding residual values for x, after 
subtracting x PM , are also depicted. 

If we adopt the assumption that GR is the correct theory of gravity for the 
PSR J 1738+0333 system, then we can get constraints on the longitude of the ascending 
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Figure 4. Probability distributions of a 2 from PSR J1012+5307 (blue dashed 
histogram), PSR J1738+0333 (red dotted histogram), and their combination (black 
solid histogram). At 95% confidence level, |&2| is constrained to be less than 1.8 x 10~ 4 
from the combined probability distribution. 

node Q. At a 95% confidence level one finds 

n E (0°, 60°) or (190°, 360°) when i < 90° , 

Q E (10°, 240°) when i> 90°. ^ ' 

We will partly use these results in the next section, where we constrain Naturally, 
to constrain a 2 we cannot use (132]) as it is based on GR, i.e. a 2 = 0. The same is true 
for PSR 1012+5307, where the constraint of Q can be found in (11-14) of }34j . 

In the calculation, for every Q 10 5 Monte Carlo simulations are implemented to 
account for the measurement uncertainties of fi a , fi$, d, v r , q, m c , and x. The results are 
plotted in the lower panels of figure |3] for i < 90° (left) and i > 90° (right), respectively 
Similarly, divergencies are caused by the unfavorable configurations with ip ~ 90° or 
~ 90° (see upper panels for reference), where a.2 can hardly be constrained. For most 
Q, a 2 is constrained to be of order ~ 5 x 10~ 5 , about two times better than that of 
PSR J1012+5307. 

3.3. Probability distribution of a.2 

We cannot directly constrain a.2 from (1251) , and have to set up Monte Carlo simulations to 
account for our ignorance in Q and ambiguity between % and 180° — i. In the simulation, 
we assume that Q is uniformly distributed between 0° and 360°, and an equal probability 
to have % < 90° or % > 90°. The measurement uncertainties of the spatial binary motions, 
the time derivative of the projected semi-major axis x, the binary mass ratio q, the 
companion mass m c , and also the distance from parallax measurement, are properly 
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accounted for, assuming a Gaussian error distribution. One million simulations for each 
of the two pulsars have been conducted. In every realization, we subtract the secular 
effects on x from proper motions, and calculate a 2 according to ( 1251) . 

The simulated results are summarized as probability distributions of ct 2 for 
PSRs J1012+5307 and J1738+0333, illustrated in figure H] using a blue dashed histogram 
and a red dotted histogram, respectively. The long tails of the probability distributions 
from these two binaries are due to the undesirable configurations when if; ~ 90° or 
$ ~ 90°, which cause the divergencies in the lower panels of figures [2] and [3j From these 
probability distributions we find at 95% confidence level, 

\a 2 \ < 3.6 x 1(T 4 for PSR J1012+5307, (33) 
\a 2 \ < 2.9 x 10~ 4 for PSR J1738+0333. (34) 

In figure HI we also show the probability distribution of a> 2 from the combination 
of these two pulsar binaries (black solid histogram), assuming that their measurements 
are independent, and that ct 2 has only a weak functional dependence on the NS mass in 
the range of 1.3 - 2.0 M Q . The combined probability distribution demonstrates a much 
shorter and suppressed tail, which means it is very unlikely that both systems are in 
the unfavorable configurations. From the combined probability distribution, we obtain 
a constraint of 

\a 2 \ < 1.8 x 10" 4 (95% C.L.) . (35) 

It is by three orders of magnitude better than the result of [51], i.e. (HI), although one 
has to keep in mind that the double pulsar tests the gravitational interaction of two 
strongly self-gravitating objects. This limit is still by three orders of magnitude weaker 
than the Solar system limit ([2]), but accounts for possible strong-field deviations in NSs. 
Compared with the Solar system limit from LLR [39], (1351) is still 3.6 times weaker. 
Being related to a secular effect, the limit on a 2 will improve fast with observing time 
T bs, namely proportional to T ob ^ 2 . A disadvantage of this test, as compared for instance 
to the LLR experiment, is its dependence on probabilistic considerations with respect 
to the unknown angle Q. 



4. A robust method to constrain c\\ 

As mentioned in section [21 Damour and Esposito-Farese used a novel geometrical way 
to constrain 6l\ with small-eccentricity binary pulsars p2]. In their paper they showed 
that the observed eccentricity vector, e(t), is a vectorial superposition of a "rotating 
eccentricity" e#(i) (with constant length), and a "forced eccentricity" ep, 

e(t) = e F + e R (t) , (36) 

where e^(t) is a vector of (unknown) constant magnitude which rotates in the orbital 
plane with angular velocity cOpn, and is a fixed vector, 
«i q ~ 1 n b 

e F = y^> —T - — V k x w . (37) 
4<r q + 1 a; PN 
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Figure 5. Illustration of our robust method to constrain &i. Upper: (i) observational 
eccentricity e(i), is a vectorial superposition of a "rotating eccentricity" e#(i) and a 
"forced eccentricity" ep [12]; (ii) the worst starting configuration, where e#(i) rotates 
A# = 26 during the observational time, and in the midtime, en(t) and ep cancel 
out completely; this applies to binaries where no explicit eccentricity is detected, 
like in PSR J1012+5307, and it constrains &i most conservatively; (iii) during the 
observational time, en(t) rotates out an angle A9, and the time-averaged e and its 
variance should be consistent with measurements (see text for two criteria) ; this applies 
to binaries where an explicit eccentricity vector is detected, like in PSR J1738+0333. 
Lower: (a) the sinusoidal evolution of e 2 (t) as a function of 6(t) = 9q + uipj<st; the 
indicated A8 corresponds to the one of the upper panel (ii); (b) the (not sinusoidal) 
time evolution of e(t) as a function of 0(t), and the indicated A8 corresponds to the 
one of the upper panel (iii). 

A graphical illustration of this dynamics is given in the upper panel (i) of figure EJ 

In analogy with the equivalence principle violation test in a probabilistic 

reasoning is used in [12] to constrain a,\. Because of two unknown angles, i.e., the time- 
dependent angle 9{t) between ep and e#(t), and the unknown angle of the longitude of 
the ascending node, Q, the probabilistic method has to assume a uniform distribution 
of 6(t) and Q between 0° and 360°. To make the "random #(£)" argument plausible, the 
binaries used in the test should be old enough to let e R (t) have rotated several cycles 
[T2T ITT] . Above two restrictions can be dropped in our robust method below. 

In order to address a potential selection effect in the test of [12], when picking 
the system with the best figure of merit from a whole population of binary pulsars, an 
ati test has been developed in [50] that extends the method of [12] to the full related 
population, including those systems that have a low figure of merit. However, from the 
viewpoint of alternative gravity theories, the &i parameter might depend on the specific 
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masses of the binary components used in the analysis. Such a dependence is well known 
for the generalized PPN parameters 7 and $ [52]. But in an analysis that combines 
different binary pulsar systems, one is forced to a priori assume a weak dependence of 
&i on the NS mass. Our new robust method can also naturally overcome the issues 
concerning the selection effect and the mass dependence, as it is based on a direct test 
of secular changes caused by &i in individual systems. 

The basic idea behind our robust method is the expected change in the eccentricity 
vector during the observational span T Q b s , in case of a non- vanishing a±. Because of the 
considerable periastron advance (see tabled] for GR values), and the fact that, by now, 
the timing observations span more than 10 years for PSRs J1012+5307 and J1738+0333, 
the "rotating eccentricity" e#(i) has already swept out a sizable angle A6 (> 10° for 
both PSRs J1012+5307 and J1738+0333). In the presence of a large ep, this would 
induce an observable change of the eccentricity vector (see figure |5]). We can use the 
tight constraints on the very small eccentricities of these two systems to limit such a 
variation of the eccentricity vector, and use this to directly constrain 6t\ with remarkable 
precision. 

We will distinguish between two cases, depending on whether an upper bound or 
a positive measurement of the eccentricity vector is made from observations. If the 
eccentricity is not measured (like for PSR J1012+5307), we use the measured (small) 
limit on the eccentricity to constrain ep, and consequently 6l\. If a positive detection 
of the eccentricity vector is made (like for PSR J1738+0333), we include its directional 
information along with the smallness of its variation to constrain a±. 

4.1. PSR J1012+5307, a short orbital period system with an unmeasured eccentricity 

For PSR J1012+5307, at present there is no measurement of an orbital eccentricity, but 
a tight upper limit of the order of 10 (see table [1]). This low limit, in combination 
with the fact that the periastron should have precessed by about 10° over the observing 
time span T Q b s , allows to put constraints on an c^-related polarization of the orbit. In 
the following we will outline the method. 
From (136|) . we arrive at, 

e 2 = e 2 F + e 2 R + IcfCr cos9 . (38) 

Since cr and cf are both constant, the observational e 2 [#(£)] changes as a sinusoidal 
function of 9(t), as shown in the lower panel (a) of figure [5j As 6{t) = 9 + cj PN t is a 
linear function of time t, e 2 (t) is a sinusoidal function of t as well. We also show a typical 
temporal evolution of e [#(£)] in the lower panel (b) of figure [51 which is not sinusoidal. 

For relativistic small-eccentricity binary systems, especially for those that are being 
observed for more then a decade, the vector e#(t) has already swept out a non- negligible 
angle, which corresponds in the lower panels of figure [5] to a span A9 in the horizontal 
axis. This factor would induce an eccentricity variance because of the vectorial addition 
in ( 136|) that scales with ct\. If then over a long time span an orbital eccentricity remains 
undetected to a small value, like for PSR J1012+5307, one can directly constrain 
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to some upper limit |ai| upper , without any probabilistic assumptions about 9. The 
reasoning is as follows. If ep 3> then e ~ ep, hence the observational smallness of e 
directly puts a limit on e^. If <C e^, then e ~ e#, hence ~ e. Therefore, 

a comparable magnitude of and e# which cancels them out (for an appropriate 9) 
would present the weakest constraint on e^. But even in the unlikely event of complete 
cancellation, a sizable di cannot hide forever. A finite observational time T Q ^ S , hence a 
finite change in 9, A9 = cupN^obs; would induce a sizable variation of the eccentricity. 
We can use this "induced" eccentricity variation to perform the most conservative limit 
on oli. 

As we can see in the lower panel (a) of figure El for a given observational time span, 
the most conservative configuration is the one in which ep and e^(t) cancel out right 
at the middle of the observational time span, as illustrated in the upper panel (ii). For 
the most conservative configuration, the "rotating eccentricity" sweeps out an angle 9 
from Ti — 5 to n + 5 during the observational time span T obs , where A9 = up^T ohs = 25. 

Before moving on, we would like to point out that our analysis accounts for the 
fact that the timing eccentricity published results from a fit to the whole observational 
data set of the pulsar binary, spanning 15 years. Hence the published eccentricity (table 
[1}) represents a "weighted" average of a potentially changing eccentricity e(t). After 
accounting for this, we can get an upper limit for the maximum eccentricity e hidden 
in the data, and from this a limit on ej? from ( 1381) . 



ep - h 6 ffm ' (39) 

a/1 - sm z S/S 2 

which can be converted into an upper limit on di through (l37l) . 

la.pp- = -^JL^ !3? ■ *-"„ = ( *») (S.Y 1 , (40) 



vr 2 (g-l) x y/i - sin 2 5/5 2 \ sin ^/ 

where ip is again the angle between the orbital angular momentum and w (see figured]). 
On the right hand side of (140jl . for a given Q, all quantities are observables or can be 
directly derived from observational quantities, except up^. As given in f JTOj) . cj PN is the 
advance rate of periastron, with potential corrections from the (generalized) Eddington- 
Robertson-Schiff parameters and PFE parameters. The solution is to take advantages of 
the smallness of 5, as it is indeed the case in the binaries which we will use to constrain 
&i. For PSRs J1012+5307 (and J1738+0333) the 5 is smaller than 10°. When 5 is 
small, PU|) becomes, 

I upper _ 2^ eP b 2 ( sini \ (w 



7T 2 (q - 1) xT ohs \smipj \cJ 
where, to first order, wpn cancels out in the numerator and denominator, since 
y/l — sin 2 5/5 2 ~ 5/\/3 = ^T obs /(2v / 3). Strictly speaking, the smallness of 5 for 
PSR J1012+5307 (and PSR J1738+0333) has been inferred from the GR value of w PN . 
The argument would break down if there is a factor of a few deviation from GR in these 
systems. This we consider as unlikely, as there is neither such a deviation from GR in 
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Figure 6. Upper limit of |di| as a function of the longitude of ascending node f2, 
derived from PSR J1012+5307 by using flU]) (Upper: i < 90°; Lower, i > 90°). The 
shadowed regions are the allowed values of Q at 95% confidence level, assuming 612 = 
[34] . Different contours correspond to 68%, 95%, and 99.7% confidence levels. 



the gravitational wave damping of these systems, nor is such a large deviation seen in 
generic direct tests of oj, like in the double pulsar [32] . 

As mentioned in section [3j PSR J1012+5307 has been observed for 15 years (for 
the timing solution presented in [M]). During that time, its periastron has already 
rotated out an angle A6 ~ 10°, which corresponds to 5 ~ 5° in the most conservative 
configuration. The smallness assumption of S is satisfied, hence we can use ( I4"TT) instead 
of the more rigorous expression (140|) . 

After properly accounting for all measurement errors, we carry out 10 5 Monte Carlo 
simulations to get an upper limit of a.\ for every value of Q between 0° and 360° (in 
steps of one degree). The results are depicted in figure [6] for i < 90° (upper) and i > 90° 
(lower), respectively. Our calculation uses the worst configuration (see the upper panel 
(ii) of figure |5]), hence the limit is most conservative and reliable. It is easily seen that 
the results for i < 90° and i > 90° are merely shifted by 180°. The existence of the 
peak near 15° for the case i < 90° (195° for the case % > 90°) is caused by the 1/ sinijj 
factor in ( 140]) . It can be understood from the ip curve in the upper panels of figure |2j 
It is important that sin^ does not vanish for Q G [0°,360°), as this avoids divergencies 
like in figure |2] for the a.2 test, which would have to be excluded based on probabilistic 
considerations. 

From figurelHl most Q realization would limit |&i| to be less than 2x 10~ 4 (95% C.L.). 
Worth to mention that, our confidence level for the &i test is purely from measurement 
errors, in contrast with that from probabilistic assumptions. If we assume a random 
fl G [0°,360°), a similar constraint is obtained. However, since we want to constrain a\ 
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most robustly, we conservatively adopt the worst configuration (Q ~ 15° for % < 90° or 
ft ~ f95° for i > 90°) and get a limit 

|«i| < 1.3 x 1(T 3 (95% C.L.) . (42) 

This limit is one order of magnitude worse than that obtained in [50], see ([3]). But as it 
avoids the probabilistic considerations of the method used in [50], we consider this limit 
as more robust. Furthermore, this limit is more likely to improve in the future than 
that of [50], which we will discuss when giving the figure of merit of this test, in section 

It is worth mentioning that, as we can see in section [21 for small-eccentricity binaries, 
the effects induced by &i and a 2 decouple. Hence this kind of test is not directly 
influenced by a non-zero &2- But a vanishing &2 would tighten the &i constraint a little 
further. If we adopt a zero &2, or take the Solar limit ([2]) for a 2 as a limit for & 2 , then 
the observed x of [31] can be attributed totally to the contribution of the proper motion, 
i.e. (I29I) . Consequently, as mentioned, Q can be constrained to certain value ranges [34|. 
We plot the consistent Q values (95% C.L.) as shadowed regions in figure El We can see 
that, this extra constraint excludes the worst Q configuration in both cases of i < 90° 
and % > 90°. Hence, with a vanishing d 2) we get a much tighter (conservative) limit of, 

|di| < 1.6 x 10" 4 (95% C.L.) . (43) 

This limit is comparable to the current best Solar system limit on a\, coming from LLR 
[39], and is only slightly worse than the current best limit ([3]) for strongly self-gravitating 
bodies, but based on a method that avoids probabilistic considerations in terms of the 
exclusion of certain unfavorable angles. 

4-2. PSR J1738+0333, a short orbital period system with a measured eccentricity 

Like PSR J1012+5307, PSR J1738+0333 is also a short orbital period NS-WD binary 
which can be use to constrain ci\. This binary has been observed for 10 years [20], and 
during this time, the periastron has advanced by A8 ~ 16°. Unlike PSR J1012+5307, 
PSR J1738+0333 has a 3-a measurement of the (intrinsic) orbital eccentricity vector 
(see rj and k in table []]), while any change in the eccentricity vector is still hidden 
in the measurement uncertainties. In a PFE test that exploits all the available 
information about PSR J1738+0333, this fact has to be taken into account (cf. [48J, 
where such a directional information, i.e. the longitude of periastron u, has been used 
to improve the constraints on a violation of the strong equivalence principle). But 
let us first, for simplicity, make only use of the upper limit for the eccentricity (2- 
o~ upper limit for e: 5.7 x 10~ 7 ) and apply the method discussed in the previous 
subsection on PSR J1012+5307, as this method is fast and can be easily compared to 
the PSR J1012+5307 result, which is based on the same method. For PSR J1738+0333 
one finds S ~ 8°, and furthermore 

|«i| < 1.6 x 10~ 4 (95% C.L.) , (44) 
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Figure 7. The &\ constraint from PSR J1738+0333, as a function of the longitude 
of ascending node ft (Upper: i < 90°; Lower, i > 90°). The limit is obtained from 
Monte Carlo simulations, by using full information of the observed eccentricity vector. 
The shadowed regions are the allowed values of f2 at 95% confidence level, assuming 
a 2 = 0, see (|32|). Different contours correspond to 68%, 95%, and 99.7% confidence 
levels. 



for the worst configuration of Q and 6* (i.e. upper panel (ii) of figure[5]) , without assuming 
a vanishing a^. This limit has to be compared to ( l42i) . The improvement over that of 
PSR J1012+5307 comes from the shorter orbital period and a smaller 1/ sin^ factor for 
the worst fl, as can be seen in the upper panels of figure [3j 

However, for PSR Jl 738+0333 we have a 3a-measurement for the eccentricity 
vector, instead of just an upper limit on its magnitude (see table [p. The result is 
expressed in terms of the first and second Laplace-Lagrange parameters, rj = esino; 
and k = e cos a; [20], from the ELL1 timing model [33]. Therefore, essentially we have 
constraints on the magnitude and the direction of the eccentricity vector. To fully 
include this information, we develop a new method, whose basic idea is depicted in the 
upper panel (iii) and the lower panel (b) of figure El and is to be elaborated below. 

We a priori have no knowledge about the magnitudes of ep and e#, nor the initial 
angle between them, # , at the time when timing observation started. Instead, we have 
rough information about the superposed eccentricity e, including its magnitude and 
direction, to ~ 3cr precision (see r] and n of PSR J1738+0333 in table 1), and also the 
direction of ep as a function of the unknown Q. However, we cannot use e directly, 
for rj and n were treated as constants when fitting to TOAs, therefore they represent 
equivalently time-averaged quantities of the observational span. We set up Monte Carlo 
simulations to select those a\ which do not conflict with the timing observation of PSR 
J1738+0333. Our simulation mainly includes the picking of an &i (hence e^) and e# (a 
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sufficiently large range for > and a 9 G [0°,360°)), for every Q £ [0°,360°). The 
worst constraint of &i is announced as our most conservative limit. All measurement 
errors are properly considered in our simulations. 

The practical implementation to get the constraints for &i is as following. First, 
we randomly choose &i, Cr (> 0), uniformly from wide ranges, and 9q from [0°, 360°), in 
order to map out a data-cube for these three values. Second, for every point in the data 
cube we let the eccentricity vector, e(t) = ep + e^(t), evolve under the dynamics given 
by the Lagrangian for a time span of observations. Then we calculate the time-averaged 
7] and k, as well as their variances, a v and o~ K . We impose the following two criteria 
to get possible &\, en, 9q combinations which can still hide in the timing data: 1) The 
time-averaged r\ and k should agree with observational values within 2a measurement 
errors (we check that the results are practically unchanged if we use la or 3a instead); 2) 
The intrinsic variances of these two parameters, induced by their time evolution, should 
be smaller than the observational errors, because otherwise, they would contradict the 
actual measurement of 77 and k. By this, we assume that the reported measurement 
errors are the squared addition of the intrinsic variances and other possible errors, 
including those possibly from measurement devices and timing models. For each f2, we 
accumulate 10 5 events in total, and get a distribution from it. The median values and 
distribution widths of di are reported in figure [7]for i < 90° (upper) and i > 90° (lower), 
clS 8b function of f2. 

We can see from figure [7] that, a± is constrained to the level of ~ 10 -5 . More 
importantly, in contrast to PSR J1012+5307, the limit only weakly depends on the 
(presently) unknown angle of the ascending node, Q. For the worst configurations 
(Q ~ 92° for i < 90° and SI ~ 273° for i > 90°) one finds 



which is about 40 times stronger than the limit from PSR J1012+5307 (cf. (j4"2"|) ). As 
outlined above, this limit is free of any probabilistic considerations related to unknown 
angles. The limit in (1451) is more than five times better than the present best limit on 
«i, coming from LLR [39]. Furthermore, it is also about four times better than the less 
robust test of [50] . 

Like in the case of PSR J1012+5307, we could adopt the very small Solar system 
limit for a 2 as a limit for 0,2 , in order to constrain the range of Q with the help of x. 
But here this would only slightly improve compared to the above constraint, because of 
the weak dependence on Q. 

It is also instructive to extract from our calculations the possible values of ej?, in 
comparison to the observed eccentricity e. We find in our Monte Carlo simulations that 
ep < 1.4 x 10 -6 (95% C.L.) for the most conservative configuration in terms of i and Q. 
Compared to the (Shapiro corrected) observed e = \J"i] 2 + k 2 = 3.4 x 10~ 7 , this is only 
a factor of a few larger. This fact reveals a generic feature of our test, that a too large 
ep cannot hide in the timing data. It lays down the reasoning of the tt\ test presented 
here, utilizing long-term pulsar timing of small-eccentricity binaries. 




(45) 
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The method developed here for PSR J1738+0333, is capable of incorporating the 
directional information of a measured eccentricity vector to constrain a.\. In principle 
it can also be applied to binary pulsars where no eccentricity is measured, and the 
direction of the eccentricity vector is unconstrained. For instance, for PSR J1012+5307 
this method yields 

£i = -0.0iJ;J x 1(T 3 (95% C.L.). (46) 

As expected, this limit is comparable with that of (1421) . and constitutes a nice test of 
the implementation used PSR J1738+0333. In general, for systems that do not have a 
measured eccentricity, the method outline in section 14.11 is preferable, as the method of 
this subsection is computationally considerably more expensive. 

5. Discussions and summary 

5.1. Constraints on PFEs from the Galactic frames 

When using the isotropic CMB frame as the assumed preferred frame, we are basically 
assuming that the preferred frame is determined by the global matter distribution in 
the Universe, and that the extra vectorial or tensorial components of gravitational 
interaction are long range, at least comparable to the Hubble radius. While this is 
generally the most plausible assumption, it is still interesting to consider other, more 
local preferred frames, like the one related to the rest frame of our Galaxy, or a frame in 
Galactic co-rotation with the overall local matter at the pulsar's location (cf. [H]). It is 
straightforward to apply the computations of this paper to these two Galactic frames. 

We used the Galactic model of [44], which assumes a distance of 8.0 kpc between 
the Solar system and the Galactic center, to extract the Galactic rotation curve. Then 
the velocities of the pulsar binary with respect to the Galactic frame and the local co- 
rotating Galactic frame are obtained. They are typically smaller than the velocity with 
respect to the CMB frame (see table [2]). The binary velocity with respect to the local 
co-rotating frame is about 100 km s -1 , which is, however, one order of magnitude larger 
than that of the Solar system [38]. Hence, pulsar binaries have the advantage to probe 
secular PFEs with respect to locally co-rotating Galactic frames, thanks to their peculiar 
velocity produced by the supernova, while Solar system tests are expected to be clearly 
less sensitive to such PFEs, because of the small peculiar Solar velocity^ Table [2] gives, 
besides the limits for the CMB frame, the limits for the two Galactic frames introduced 
above, based on the methods presented in sections [3] and HI The limits concerning the 
two Galactic frames are somewhat weaker than that of the CMB frame, because of the 
smaller peculiar velocities of the pulsar binaries. Nevertheless, they could be of interest 
for tests of specific gravity theories that might predict or even require (to pass Solar 
system tests) vector or tensor fields which are aligned with the Galactic or local matter 
distribution. 

% Nevertheless, see [47] for a constraint of ot\ = (1.6 ±8.0) x 10~ 3 (95% C.L.) from LLR when choosing 
the Barycentric Celestial Reference System as the preferred frame, based on periodic effects. 
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Figure 8. Sky coverage in the di test by PSRs J1012+5307 and J1738+0333. The blue 
curves mark the directions in the sky where PSR J1012+5307 is insensitive, and the 
red curves those directions where PSR J1738+0333 is insensitive, due to the unknown 
f2. The sky is plotted in an Hammer- Aitoff projection using Galactic coordinates / and 
b. The longitude I increases from right to left, from I = —180° to I = +180°, while the 
latitude b runs from b = —90° to b = 90° from bottom to top. The grid gives steps of 
30°. The plot is based on a w = 369 km s^ 1 for the velocity of the SSB with respect 
to a potential preferred frame. The label "CMB" denotes the direction of motion with 
respect to the CMB frame. 



Finally, it is interesting to see the sky coverage of the two binary pulsars in terms of 
a sensitivity towards a preferred frame. In the spirit of a "PFE pulsar antenna array" , 
proposed in [51], two suitable binary pulsars can probe for a preferred reference frame 
in (almost) any direction in the sky. Figure [8] shows the combined sky coverage in the 
ai test for PSRs J1012+5307 and J1738+0333, for w = 369 km s" 1 , the same magnitude 
as the velocity of the SSB with respect to the CMB. The figure nicely illustrates why 
PSR J1738+0333 is more suitable to test the CMB frame, since the "CMB" dot is far 
away from any of the red curves, meaning that the CMB direction lies in the area where 
the system is particularly sensitive. The absence of any divergencies in figure [7] is a 
consequence of this. 



5.2. Strong-field modifications 

When comparing the results of this paper, obtained from pulsar- WD systems, with Solar 
system experiments, one has to keep in mind that alternative gravitational theories, in 
general, are expected to predict strong-field modifications of the PPN parameters due 
to the strong internal gravitational field of the pulsar. As an example, in scalar-tensor 
gravity the PPN parameter 7 generalizes to 
* _ , 2a A a B 

7 = Iab = 1 - 7- , (47) 

1 + a A OL B 

for a binary pulsar system, where a a and a B are the effective scalar coupling constants of 
pulsar and companion, respectively [T3]. The weak- field PPN parameter 7 is recovered 
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for a a = 01 b = OiQ. In GR one has 7 = 7 = 1. In the strong-field regime of a NS, 7 
can deviate significantly from 7 due to strong- field scalarization effects [Hj. Similarly, 
we may expect that 6l\ and oti deviate from their PPN correspondents, a\ and a^- In 
the absence of non-perturbative effects, one can illustrate this as an expansion in the 
compactnesses ca and cb of the bodies [15]. In our case, we would write something like 



where a = 1,2 and /Q and /Qj are coefficients characterizing deviations from general 
relativity, and q ~ Grrii/RiC 2 with mass rrii and radius Ri of body i. The compactnesses 
for the Earth and the Sun are roughly c e ~ 10~ 10 and c ~ 10~ 6 , respectively, 
which suppress /Q- and /Q.,- related physical effects dramatically. In contrast, NSs have 
cns ~ 0.2, which is one of the reasons why pulsar timing experiments are ideal probes 
for gravity effects associated with strong gravitational fields. Consequently, in a NS-WD 
system (cwd ~ 10~ 4 ) we could still have a significant a a , even if there is a tight Solar 
system constraint for a a . This, for instance, supports the importance of the 612 limits 
obtained in this paper. 

Finally, when discussing the constraints on parameters of alternative gravity 
theories, one should be aware of a potential compactness-dependent (or mass- 
dependent) nature of these parameters. Especially when combining different systems, 
like the c\\ test in [50], or the combined probability distribution function of the 612 test 
presented in this paper (black solid histogram of figure Hj). Such tests implicitly assume 
that the parameter is approximately the same for all systems under investigation. We are 
aware of this potential problem in our calculation. In the case of the PSRs J1012+5307 
and J1738+0333 binary systems one can argue that the similarity in the masses justifies 
such an assumption. However, in the presence of phenomena related to some critical 
mass, like the spontaneous scalarization discovered in [H], even a small difference in 
masses does not allow such an assumption. Our proposed robust test for &i overcomes 
this "mass-dependence" problem by only considering one system. Hence our final results 
on cti are more suitable to be quoted along with mentioning the specific system were it 
was obtained from and the related neutron-star mass. 

5.3. Figures of merit and further potential improvements 

Finally, let us discuss potential improvements of the current limits of a\ and &2- For 
this, one conveniently identifies the figures of merit for the different tests. Besides 
details related to the shape and the size of the orbit, various geometrical angles and 
the barycentric velocity with respect to the preferred frame (w) play a role in our 
tests. Therefore, in principle our figures of merit would depend on the geometrical 
configuration of binary systems under consideration, as well as their proper motions with 
respect to the preferred frame. However, after we drop the geometrical dependencies 
and ignore potential difference in w, we can roughly get a figure of merit of our tests. 

As for the &i test, the traditional method has a figure of merit, P^ 3 ' je [12] . 
which means that the strength of this test is not to improve until new systems with 



Ota OL a ~\~ K^iCi -\~ K^ijCiCj -]-..., 
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higher P b je are discovered. In contrast, our figure of merit for the robust 6l\ test is 
^obs/(-Pfe 4//3 e) for binary pulsars with unmeasured eccentricities e < e (cf. (14TT)). Besides 
the discoveries of new systems with smaller P^ 3 e, the constraint has the potential 
to improve when the observational span becomes longer. In fact, it improves as 
T^ s 2 , as long as e remains smaller than e, and ooT ohs < 1. For a binary pulsar with 
measured eccentricity e, we get a similar figure, where e is to be replaced with a e , the 
measurement error of the eccentricity vector, or that of the first and second Laplace- 
Lagrange parameters, namely a v and a n . If in the future the secular evolution of the 
eccentricity vector in PSR Jl 738+0333, due to the relativistic periastron advance, can 
be measured, we could further constrain a potential polarization of the orbit caused by 
a non-vanishing a±, or even detect the presence of a significant PFE eccentricity ep. 

In contrast to our secular effects, the limit from LLR is based on periodic effects 
(see e.g. p2]), and therefore only improves as T Q ^ S . Moreover, to test PFEs in LLR 
one has to deal with the motion of the Earth-Moon system around the Sun, introducing 
tidal forces and an annually changing w, hence leading to the theoretical complexity of 
a three-body problem in the presence of a preferred frame |17j . 

For the a 2 test one finds the figure of merit to be l/(P^ 3 a±) from ( |25|) . where 
a± is the measurement uncertainty of x. Hence, more relativistic systems (smaller Pb) 
with high timing precision (especially of x) are advantageous to do the ct^ test. For the 
systems discussed here, namely PSRs J1012+5307 and J1738+0333, persistent timing 
observations will reduce the measurement uncertainty of o~±. Hence the &2 test will 
improve continuously as T 3 ^, in contrast with that of [12] ■ An even faster improvement 
is expected to come from new receiver and backend technologies and new telescopes. 

It is worth noting that a measurement of the unknown longitude of the ascending 
node, Q, would improve both the 6t\ and the a<i test. If Q can be determined 
independently, even rough constraints on Q would make these tests more efficient. 
It would eliminate the (systematical) "double peak" structure in the probability 
distribution function of &2 in figure HI Also it would select a specific limit on 6t\ from 
figures [6] or [71 where presently we are conservatively using the worst Q configuration. 
Unfortunately, neither pulsar timing nor optical astrometry are likely to provide such 
a measurement in the near future. Maybe scintillation measurements would be able to 
provide interesting constraints on fl, like this is the case for the double pulsar (Rickett 
et al, in prep.). 

5.4- A brief summary of the results 

In summary, we presented an extended orbital dynamics of pulsar binaries under the 
influence of preferred frame effects that accounts for both generalized PPN parameters, 
&i and «2- In the limit of a small eccentricity, orbital effects from a\ and decouple. 
We implemented two new methods to constrain 6t\ and «2, by directly constraining 
secular orbital changes expected from a violation of local Lorentz invariance in the 
gravitational sector. Both methods have been applied to the two binary pulsars PSRs 
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J1012+5307 and J1738+0333, where we have full spatial velocity information. For a 
frame at rest with respect to the CMB, the best limit we obtain is 

a t = -OAtli x 10" 5 (95% C.L.) , (49) 

which avoids the probabilistic considerations of previous methods, and clearly surpasses 
the current best limits obtained with both, weakly (Solar system) and strongly (binary 
pulsars) self-gravitating bodies, namely (CD) and ([3]). 
Concerning a 2 , the best limit we obtain is 

\a 2 \ < 1.8 x 10" 4 (95% C.L.) . (50) 

This limit is still three orders of magnitude weaker than the weak-field limit obtained 
in the Solar system, but constrains possible deviations related to the strong internal 
gravitational fields of NSs. The limit here surpasses the current best limit for strongly 
self-gravitating bodies, namely (T4|), by three orders of magnitude, although strictly 
speaking they are different in their physical nature, as (jlj) probes the interaction between 
two strongly self-gravitating bodies, in contrast to the pulsar- WD systems used in this 
paper. A drawback of the a 2 limit presented here is that it is still based on probabilistic 
considerations in excluding unfavorable values of the longitude of the ascending node 
Q, and the combination of two systems with different NS masses. 

Our new methods promise continuous improvements with on-going timing 
observations of known systems, as well as the discovery of new suitable systems. In 
addition, new receiver and backend technologies as well as new telescopes, like FAST 
[4*0] and SKA [45], will tremendously improve these tests. 

Acknowledgment s 

We thank John Antoniadis, William Coles, Paulo Freire, Michael Kramer and Clifford 
Will for stimulating discussions. We are grateful to Michael Kramer for carefully reading 
the manuscript. Lijing Shao is supported by China Scholarship Council (CSC). 

References 

[1] J. Antoniadis, M. H. van Kerkwijk, D. Kocster, P. C. C. Freire, N. Wex, T. M. Tauris, 
M. Kramer, and C. G. Bassa. The relativistic pulsar-white dwarf binary PSR J1738+0333 
I. Mass determination and evolutionary history. MNRAS, 423:3316, 2012. 

[2] Z. Arzoumanian, K. Joshi, F. A. Rasio, and S. E. Thorsett. Orbital parameters of the PSR 
B1620— 26 triple system. In S. Johnston, M. A. Walker, & M. Bailes, editor, IAU Colloquium 
160 - Pulsars: Problems and Progress, Astronomical Society of the Pacific Conference Series, 
105:525, 1996. 

[3] Q. G. Bailey and V. A. Kostelecky. Signals for Lorentz violation in post-Newtonian gravity. 

Phys. Rev. D, 74:045001, 2006. 
[4] B. M. Barker and R. F. O'Conncll. Relativistic effects in the binary pulsar PSR 1913+16. ApJL, 

199:L25, 1975. 

[5] J. D. Bekenstein. Relativistic gravitation theory for the modified Newtonian dynamics paradigm. 
Phys. Rev. D, 70:083509, 2004 [Erratum: Phys. Rev. D, 71:069901, 2005]. 



New tests of LLI of gravity with binary pulsars 



30 



[6] J. F. Bell, F. Camilo, and T. Damour. A tighter test of the local Lorentz invariance of gravity 

using PSR J2317+1439. ApJ, 464:857, 1996. 
[7] M. Burgay, N. D'Amico, A. Posscnti, R. N. Manchester, A. G. Lyne, B. C. Joshi, 

M. A. McLaughlin, M. Kramer, J. M. Sarkissian, F. Camilo, V. Kalogera, C. Kim, and 

D. R. Lorimer. An increased estimate of the merger rate of double neutron stars from 

observations of a highly relativistic system. Nature, 426:531, 2003. 
[8] P. J. Callanan, P. M. Garnavich, and D. Koestcr. The mass of the neutron star in the binary 

millisecond pulsar PSR J1012+5307. MNRAS, 298:207, 1998. 
[9] S. Chandrasekhar and G. Contopoulos. On a post-Galilean transformation appropriate to the 

post-Newtonian theory of Einstein, Infold and Hoffmann. Royal Society of London Proceedings 

Series A, 298:123, 1967. 

[10] T. Damour and N. Deruelle. General relativistic celestial mechanics of binary systems. II. The 

post-Newtonian timing formula. Ann. Inst. Henri Poincare Phys. Theor., 44:263, 1986. 
[11] T. Damour and G. Schafer. New tests of the strong equivalence principle using binary-pulsar data. 

Physical Review Letters, 66:2549, 1991. 
[12] T. Damour and G. Esposito-Farese. Testing local Lorentz invariance of gravity with binary-pulsar 

data. Phys. Rev. D, 46:4128, 1992. 
[13] T. Damour and G. Esposito-Farese. Tensor-multi-scalar theories of gravitation. Class. Quantum 

Grav., 9:2093, 1992. 

[14] T. Damour and G. Esposito-Farese. Nonperturbative strong-field effects in tensor-scalar theories 
of gravitation. Physical Review Letters, 70:2220, 1993. 

[15] T. Damour and G. Esposito-Farese. Testing gravity to second post-Newtonian order: A field- 
theory approach. Phys. Rev. D, 53:5541, 1996. 

[16] T. Damour and J. H. Taylor. On the orbital period change of the binary pulsar PSR 1913+16. 
ApJ, 366:501, 1991. 

[17] T. Damour and D. Vokrouhlicky. Testing for gravitationally preferred directions using the lunar 

orbit. Phys. Rev. D, 53:6740, 1996. 
[18] T. Damour and J. H. Taylor. Strong-field tests of relativistic gravity and binary pulsars. 

Phys. Rev. D, 45:1840, 1992. 
[19] B. Z. Foster. Strong field effects on binary systems in Einstein-aether theory. Phys. Rev. D, 

76:084033, 2007. 

[20] P. C. C. Freire, N. Wex, G. Esposito-Farese, J. P. W. Verbiest, M. Bailes, B. A. Jacoby, M. Kramer, 
I. H. Stairs, J. Antoniadis, and G. H. Janssen. The relativistic pulsar-white dwarf binary PSR 
J1738+0333 II. The most stringent test of scalar-tensor gravity. MNRAS, 423:3328, 2012. 

[21] G. Hinshaw, J. L. Weiland, R. S. Hill, and et al. Five-year Wilkinson Microwave Anisotropy Probe 
observations: data processing, sky maps, and basic results. ApJS, 180:225, 2009. 

[22] T. Jacobson and D. Mattingly. Gravity with a dynamical preferred frame. Phys. Rev. D, 
64:024028, 2001. 

[23] B. A. Jacoby. Recycled pulsars. PhD thesis, California Institute of Technology, California, USA, 
2005. 

[24] N. Jarosik, C. L. Bennett, J. Dunkley, and et al. Seven-year Wilkinson Microwave Anisotropy 
Probe (WMAP) observations: sky maps, systematic errors, and basic results. ApJS, 192:14, 
2011. 

[25] K. J. Joshi and F. A. Rasio, Distant companions and planets around millisecond pulsars. ApJ, 
488:901, 1997. 

[26] V. M. Kaspi, M. Bailes, R. N. Manchester, B. W. Stappers, and J. F. Bell. Evidence from a 
processing pulsar orbit for a neutron-star birth kick. Nature, 381:584, 1996. 

[27] V. M. Kaspi, A. G. Lyne, R. N. Manchester, F. Crawford, F. Camilo, J. F. Bell, N. D'Amico, 
I. H. Stairs, N. P. F. McKay, D. J. Morris, and A. Posscnti. Discovery of a young radio pulsar 
in a relativistic binary orbit. ApJ, 543:321, 2000. 

[28] S. M. Kopcikin On possible implications of orbital parallaxes of wide orbit binary pulsars and 



New tests of LLI of gravity with binary pulsars 



31 



their measur ability. ApJL, 439:L5, 1995. 
[29] S. M. Kopcikin. Proper motion of binary pulsars as a source of secular variations of orbital 

parameters. ApJL, 467:L93, 1996. 
[30] V. A. Kostelecky and N. Russell. Data tables for Lorentz and CPT violation. Reviews of Modern 

Physics, 83:11, 2011. 

[31] M. Kramer, I. H. Stairs, R. N. Manchester, M. A. McLaughlin, A. G. Lync, R. D. Fcrdman, 
M. Burgay, D. R. Lorimcr, A. Possenti, N. D'Amico, J. M. Sarkissian, G. B. Hobbs, 
J. E. Reynolds, P. C. C. Freire, and F. Camilo. Tests of general relativity from timing the 
double pulsar. Science, 314:97, 2006. 

[32] M. Kramer and N. Wex. The double pulsar system: a unique laboratory for gravity. Class. 
Quantum Gray., 26:073001, 2009. 

[33] C. Lange, F. Camilo, N. Wex, M. Kramer, D. C. Backer, A. G. Lyne, and O. Doroshcnko. Precision 
timing measurements of PSR J1012+5307. MNRAS, 326:274, 2001. 

[34] K. Lazaridis, N. Wex, A. Jessner, M. Kramer, B. W. Stappers, G. H. Jansscn, G. Desvignes, 
M. B. Purver, I. Cognard, G. Theureau, A. G. Lyne, C. A. Jordan, and J. A. Zensus. Generic 
tests of the existence of the gravitational dipole radiation and the variation of the gravitational 
constant. MNRAS, 400:805, 2009. 

[35] D. R. Lorimer and M. Kramer. Handbook of Pulsar Astronomy. Cambridge University Press, 
2004. 

[36] D. R. Lorimcr, A. G. Lyne, L. Festin, and L. Nicastro. Birth rate of millisecond pulsars. Nature, 
376:393, 1995. 

[37] A. G. Lyne, M. Burgay, M. Kramer, A. Possenti, R. N. Manchester, F. Camilo, M. A. McLaughlin, 
D. R. Lorimer, N. D'Amico, B. C. Joshi, J. Reynolds, and P. C. C. Freire. A double-pulsar 
system: a rare laboratory for relativistic gravity and plasma physics. Science, 303:1153, 2004. 

[38] F. Mignard. Local galactic kinematics from Hipparcos proper motions. Astron. Astroph., 354:522, 
2000. 

[39] J. Miiller, J. G. Williams, and S. G. Turyshev. Lunar laser ranging contributions to relativity and 
geodesy. In H. Dittus, C. Lammerzahl, & S. G. Turyshev, editor, Lasers, Clocks and Drag-Free 
Control: Exploration of Relativistic Gravity in Space, Astrophysics and Space Science Library, 
349:457, 2008. 

[40] R. Nan, D. Li, C. Jin, Q. Wang, L. Zhu, W. Zhu, H. Zhang, Y. Yue, and L. Qian. The five-hundred 

aperture spherical radio telescope (FAST) project. Int. J. Mod. Physics D, 20:989, 2011. 
[41] L. Nicastro, A. G. Lyne, D. R. Lorimer, P. A. Harrison, M. Bailes, and B. D. Skidmore. PSR 

J1012+5307 - a 5.26-millisecond pulsar in a 14.5-hour binary system. MNRAS, 273:L68, 1995. 
[42] K. Nordtvedt. Probing gravity to the second post-Newtonian order and to one part in 10 7 using 

the spin axis of the sun. ApJ, 320:871, 1987. 
[43] K. Nordtvedt and C. M. Will. Conservation laws and preferred frames in relativistic gravity. II. 

Experimental evidence to rule out preferred-frame theories of gravity. ApJ, 177:775, 1972. 
[44] B. Paczynski. A test of the galactic origin of gamma-ray bursts. ApJ, 348:485, 1990. 
[45] R. Smits, M. Kramer, B. Stappers, D. R. Lorimcr, J. Cordes, and A. Faulkner. Pulsar searches 

and timing with the square kilometre array. Astron. Astroph., 493:1161, 2009. 
[46] I. S. Shklovskii. Possible causes of the secular increase in pulsar periods. Soviet Ast., 13:562, 

1970. 

[47] M. Soffcl, S. Klioncr, J. Miiller, and L. Biskupek. Gravitomagnetism and lunar laser ranging. 

Phys. Rev. D, 78:024033, 2008. 
[48] I. H. Stairs, A. J. Faulkner, A. G. Lync, M. Kramer, D. R. Lorimcr, M. A. McLaughlin, 

R. N. Manchester, G. B. Hobbs, F. Camilo, A. Possenti, M. Burgay, N. D'Amico, P. C. Freire, 

and P. C. Gregory. Discovery of three wide-orbit binary pulsars: implications for binary 

evolution and equivalence principles. ApJ, 632:1060, 2005. 
[49] N. Wex. A timing formula for main-sequence star binary pulsars. MNRAS, 298:67, 1998. 
[50] N. Wex. Small-eccentricity binary pulsars and relativistic gravity. In M. Kramer, N. Wex, & 



New tests of LLI of gravity with binary pulsars 



32 



R. Wielebinski, editor, IAU Colloq. 177: Pulsar Astronomy - 2000 and Beyond, Astronomical 

Society of the Pacific Conference Series, 202:113, 2000. 
[51] N. Wex and M. Kramer. A characteristic observable signature of preferred-frame effects in 

relativistic binary pulsars. MNRAS, 380:455, 2007. 
[52] C. M. Will. Theory and Experiment in Gravitational Physics. Cambridge University Press, 1993. 
[53] C. M. Will. The confrontation between general relativity and experiment. Living Rev. Relativity, 

9:3, 2006 [URL (cited on 2012-07-22): |http: / / www.livingreviews.org/rrr-2006-3 . 
[54] C. M. Will and K. Nordtvcdt. Conservation laws and preferred frames in relativistic gravity. I. 

Preferred-frame theories and an extended PPN formalism. ApJ, 177:757, 1972. 



Table 1. Relevant parameters in the PFE calculations for PSRs J1012+5307 [31133] and J1738+0333 [HI Q]. 



Pulsars 


PSR J1012+5307 


PSR J1738+0333 


Right Ascension, a (J2000) 


10 h 12 m 33?4341010(99) 


17 h 38 m 53?9658386(7) 


Declination, 5 (J2000) 


53°07'02"60070(13) 


03°33'10'.'86667(3) 


Proper motion in a, fi a (masyr -1 ) 


2.562(14) 


7.037(5) 


Proper motion in S, ^ig (masyr" 1 ) 


-25.61(2) 


5.073(12) 


Distance, d (kpc) 


0.836(80) 


1.47(10) 


Spin period, P (ms) 


5.255749014115410(15) 


5.850095859775683(5) 


Orbital period, P& (d) 


0.60467271355(3) 


0.3547907398724(13) 


Projected semi-major axis, x (lt-s) 


0.5818172(2) 


0.343429130(17) 


i] = e sinw (10 _7 ) a 


-1.4 ± 3.4 


-1.4 ± 1.1 


k = e cos a; (10 -7 ) 


0.6 ±3.1 


3.1 ±1.1 


lime derivative ol x, x (10 ss ) 


2.3(8) 


0.7(5) 


lime derivative or t% 9 r% [W ss ) 

XXdQldjl VClUClLy, (7^ IA1I1 t> J 


rn/i a \ 

5U(±4) 


1 7 C\(0 1 \ 

— 1 7 .(J(o.lJ 


AA(R\ 

QQ\0 j 




Mass ratio, q = m p /m c 


10.5(5) 


8.1(2) 


WD mass, m c (M ) 


0.16(2) 


n 1oi +0.008 
0.181_ 007 


Pulsar mass, m p (Mq) 


1.64(22) 


1 46+ - 06 


Mass function, / (M Q ) 


0.000578 


0.0003455 


Inclination, i (deg) 


52(4) or 128(4) 


32.6(10) or 147.4(10) 


f GR ) 1 

Advance of pcriastron (in GR), Wp N (degyr ) 


0.69(6) 


1.57(5) 


"Characteristic" velocity (in GR), V^ GR) (kms" 1 ) 


308(13) 


355(5) 



Intrinsic 77, after subtraction of the contribution from the Shapiro delay according to (A22) of (3*1 



Table 2. Limits on the PFE parameters di and d 2 at 95% confidence level, from the NS-WD binaries, PSRs J1012+5307 and J1738+0333 
(see text for details). 



Pulsar Binary 


Preferred Frame 


w(kms 1 ) 


OLl 




&2 




J1012+5307 


CMB 
Galaxy 

Local Galactic rotation 


477(14) 
157(15) 
114(23) 


|di| < 1.3 x 10~ 3 
di < 8.3 x lO- 3 
di < 7.5 x 10~ 3 


|d 2 
|d 2 
|d 2 


< 3.6 x 10" 

< 7.9 x 10- 

< 1.1 x 10" 


-4 
-3 
2 


J1738+0333 


CMB 
Galaxy 

Local Galactic rotation 


327(6) 
265(6) 
82(10) 


di = -OAtfl x 10" 5 
di = -0.3^o x 10~ 5 
di = +0.ltH x 10~ 4 


|d 2 
|d 2 
|d 2 


< 2.9 x 10- 

< 8.3 x 10- 

< 1.0 x 10" 


-4 
-4 
-2 


Combined 


CMB 
Galaxy 

Local Galactic rotation 






d 2 
|d 2 
d 2 


< 1.8 x 10- 

< 4.5 x 10- 

< 3.4 x 10" 


-4 
-4 
-:s 



